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Abstract 

We study the relation between the hessian matrix of the riemannian Regge action on 
a 4-simplex and linearized quantum gravity. We give an explicit formula for the hessian 
as a function of the geometry, and show that it has a single zero mode. We then use a 3d 
lattice model to show that (i) the zero mode is a remnant of the continuum diffeomorphism 
invariance, and (ii) we recover the complete free graviton propagator in the continuum 



X 



(3jT)! limit. The results help clarify the structure of the boundary state needed in the recent 

calculations of the graviton propagator in loop quantum gravity, and in particular its role 
in fixing the gauge. 

> 
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^ ! 1 Introduction 

■ Loop quantum gravity (LQG) and the spinfoam formalism [IJ are candidate background inde- 
pendent quantizations of general relativity (GR). A key open question to check the consistency 
of this approach is whether it gives the right semiclassical limit. Recently there have been im- 
portant advances towards a positive answer in cosmological models |2J, black hole physics [3], 
canonical LQG jl] and in the linearized expansion [SI E] . In particular in [6] it was studied if the 
riemannian theory admits a regime in which the low energy description of linearized quantum 
gravity emerges. The basic quantity of linearized quantum gravity is the graviton propagator. 
This quantity has been denned in [S] and the promising results of [6] show that the correct 
spacetime dependence of some of its components is indeed emerging in a well defined limit. 
The key for this mechanism is the intermediate step of linearized quantum Regge calculus. 
Linearized quantum Regge calculus has been studied in the literature [7J [8], and the correct- 
ness of its continuum limit poses the spinfoam calculations on solid grounds. Yet there are 
many open questions in the spinfoam procedure, such as extending the formalism to arbitrary 
backgrounds, dynamically computing the boundary state, and understanding the gauge fixing 
procedure. In this paper we investigate linearized quantum Regge calculus from a perspective 
that is close to the spirit of the spinfoam calculations, with the aim of shedding light on some 
of the open questions. 

Our first result is a formula for the hessian of the riemannian Regge action on an arbitrary 
n-simplex configuration. This formula generalizes the equilateral configuration case appeared in 
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the literature (see third reference in [6]), and it is relevant for studying the graviton propagator, 
or 2-point function, around a generic background. The latter scales as the background distance 
between the two points to the power 2 — n, so one might naively expect a similar scaling 
for the hessian matrix. The explicit formula we obtain shows that this is not the case, and 
furthermore shows that this matrix is not invertible. This is reminiscent of the non-invertibility 
of the quadratic term of the linearized continuum theory due to diffeomorphism invariance, and 
raises the question whether the non-invertibility of the hessian can be traced back to the same 
symmetry. 

To address both questions of scaling and non-invertibility, we consider a 3d infinite lattice 
model. In this model we show that adding up the hessian matrices of all the tetrahedra in the 
full lattice one obtains the correct quadratic term of linearized Regge calculus. This in turn 
can be directly obtained discretizing the (quadratic term of the) continuum Einstein-Hilbert 
action. Upon properly gauge fixing this quadratic term, we obtain the correct free propagator, 
including both its distance dependence and its tensorial structure. This is our second result, 
and it allows us to answer the above questions. First, the correct scaling is indeed encoded 
in the hessian matrix, and shows up explicitly when a full lattice is considered. Second, the 
non-invertibility is the remnant of the diffeomorphisms on the lattice, and the full lattice gauge 
fixing term induces a gauge fixing term for a single hessian matrix. In our analysis we benefited 
from the techniques and results of [8], which we develop further. 

Based upon these results, we then discuss in more details the gauge fixing procedure and 
the general boundary formalism. The discussion allows us to make a few important points 
concerning the spinfoam calculations of [6]. In particular, we argue that in those calculations 
the gauge is fixed by the boundary state. This means that the free parameters entering it have 
to be fixed not only by the dynamics, but also by the choice of gauge. We give a prescription 
for computing dynamically the quadratic order of the boundary state. 

Our results support once more the relevance of Regge calculus to study the semiclassical 
limit of spinfoams. This idea has been often advocated [9], and dates back to the original 
Ponzano- Regge model [10J. The link between spinfoam kernels and Regge calculus has been 
deeply studied [UJ , and recently extended in 3d to expectation values of geometry [12] and to 
the coupling with Yang- Mills theory |13j . 

This paper is organized as follows. In section [2] we briefly review conventional length Regge 
calculus in any dimension n, discuss the modification where the n — 2- volumes of the hinges are 
the fundamental variables, and provide the explicit formula for the hessian. In particular, we 
show that this matrix has a single zero mode related to the vanishing of the Gram determinant 
for a flat n-simplex. In section [3] we introduce the 3d infinite lattice model, and show how the 
hessian matrix leads to the correct free graviton propagator. In section H] we discuss in more 
details the gauge fixing procedure, and how the hessian zero mode is related to the continuum 
diffeomorphism invariance. In section [5] we discuss the general boundary formalism, and we 
stress a few key properties that the boundary state has to satisfy to provide the correct physics. 
Finally, in section [6] we summarize our results and discuss their perspective on the spinfoam 
calculations. 
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2 Hessian of the Regge action 



Regge calculus [14] is a discrete approximation to GR where the fundamental variables are the 
lengths £ e of the edges of a triangulated manifold representing spacetime, denoted A + dA 
where dA is the boundary triangulation. In any dimension n, the curvature is distributional 
and concentrated on the n — 2 simplices of the triangulation, called hinges (edges in 3d, triangles 
in 4d). The action is 

Sr[Q =^WU+We], (1) 
heA 

where Vh is the (n — 2)-volume of the hinge h and the deficit angles represent the curvature. 
They are defined as = 2n — J2 a ®hi where 9% is the dihedral angle of the n-simplex a sharing 
the hinge h. This is a function of the edge lengths of the n-simplex as described below. The 
boundary term is SQ^[l e ] = J2h Vfc J2a ( n ~ @h(le)), an d using the definition of the deficit angles, 
one can rewrite the full action as a sum over n-simplices, 

sk[q = i>*[a sm = Yl v h - mS) + E ^ (w - « ( 2 ) 

where iV^ is the number of n-simplices hinged on h. The continuum limit of ([T|) gives the 
Einstein-Hilbert action. 

The link between LQG and Regge calculus is particularly transparent in 3d: in a well defined 
limit, the spinfoam partition function for LQG is dominated by exponentials of the Regge action 
(PQ) (for details see for instance pQ). In 4d on the other hand, the same limit seems related to a 
modification of ([1]) in which the fundamental variables are the hinge areas A t [15] . Area Regge 
calculus is much less known in the literature, and the key open question is whether one can do 
the change of variables from lengths to areas in conventional Regge calculus. A positive answer 
would support area Regge calculus as an equally viable discrete approximation of GR. This 
question has been only partially addressed so far, and at the present state of investigation, two 
results are clear: 

• A single 4-simplex has ten edges and ten triangles, thus one can directly study the Ja- 
cobian matrix for the change of variables from lengths to areas, The latter turns 
out to be non singular in various regions of the configuration space [TB] (typically, it be- 
comes singular in the presence of configurations with right angles between the triangles); 
therefore in those regions the dihedral angles are well defined functions of the areas and 
it makes sense to consider the following 4-simplex action, 

S a (At) = J2 A t{ n - 9 ^)- ( 3 ) 

t 

• On a generic triangulation, the situation is more complicated: there are more triangles 
than edges, thus constraints on the area variables are needed when many 4-simplices are 
considered, 

S R (A t ) = J2Aet(A)+C(A t ). (4) 

t 
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The role of the constraints is to guarantee that a tetrahedron shared between two 4- 
simplicies has the same geometry in both sets of area variables. Regretfully, the explicit 
form of these constraints is still unknown in the literature (see however [TTl [T8l [T9l 20]). 
and pursuing the study of these constraints is an important problem that deserves future 
work. 

A key object of interest for us is the hessian matrix of To keep our discussion general, 
we consider a single simplex in any dimension n. For a single n-simplex it is convenient to 
introduce a notation based on the vertices of the n-simplex, labeled by i — 1 . . . n + 1. Then 
l e = denotes the length of the edge e linking the vertices % and j, and 9^ denotes the 
(internal) dihedral angle between the (n — l)-simplex i (i.e. the one made of all vertices but 
i) and j, with the convention 8u = ir. We also pose V to be the n-volume of the simplex, 
Vi the (n — l)-volume of i, and Vy the {n — 2)-volume of the hinge between i and jl3 As in 
4d, an n-simplex has the same number of edges and hinges, n(n 2 +1 - ) , thus the jacobian matrix 
j£f = 2pi i s well defined. In the following, we restrict our analysis to configurations where 
hinges and lengths have a non- vanishing jacobian determinant, so we can define the n-simplex 
Regge action 

^[^=^V^7T-%(10). (5) 

Using the double index notation, the dihedral angles are given by 

n - 1 Vi(l e ) Vj{l e ) 

and leiVij) is well defined being det J*- ^ 0. 

The first derivatives of (jSJ) give simply the dihedral angles, 

This crucial property can be easily proved from the n-dimensional Schlafli identity, 

Y,V ij dO ij = 0. (8) 

A direct consequence of (JTj) is that the hessian is equal to the (n(n + l)/2 x n(n + l)/2) jacobian 
matrix for changing variables from n — 2-volumes to dihedral angles, 

= d 2 s a = _ a% 

This change of variables is singular: from dSJ) we immediately see that V^j is a left and right null 
vector for (Q. This turns out to be the only null vector: below we give the explicit expression 



lr The double index notation is very useful to deal with an arbitrary dimension n of spacetime. On the other 
hand, care is required in 3d, where confusion might arise because both and Vij are variables for the edge 
lengths. Vij is the hinge between the triangles i and j, namely the edge opposite to Uj. Using the notation ij 
to denote the edge opposite to Uj, we have Vij — Ijj. This notation will be useful in section [3] below. 
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for Qijki and show that it has rank "( n+1 ) — 1 and that the vanishing of its determinant is related 
to the closure of a flat n-simplex. 

As det J® 7^ 0, we can rewrite © as 

^. = -££fe. («•) 

rnn 

This expression allows us to evaluate Q in two steps, (i) computing the matrix gj^-, and (ii) 
computing the inverse matrix of d g™ n . The latter has a simple form in any dimension (see 
appendix Rjl . being non zero only when Jjy belongs to the hinge V mn . In particular in 4d 
it has only 30 non- zero entries (see for instance [IB]), and its determinant is different from 
zero provided we exclude configurations including right angles between triangles, as mentioned 
above. Then, its inverse entering fflUj) is well defined and it can be computed algebraically. 

The first term is the most interesting part (and furthermore its analysis also concerns 
conventional length Regge calculus). It can be explicitly computed as shown in the Appendix, 

™« - tan 6 l3 A ( In V + In V v - In V, - In V 3 ) = i -±L- M C m , (11) 

where 



94/ J <%A lJ J J n 2 sin% V" 2 



CjjfcZ = ^ COS COS 9u + COS 0jjfc COS 9jlj COS % + COS Oik COS 6^7 + COS Qjk COS (12) 

Explicitly, 

(sin 2 #i2 - sin 2 9 12 cos 6> 2 3 • • • \ 

-sin 2 # 13 cos# 23 sin 2 # 13 ... . (13) 

Notice that the matrix C^i is not symmetric under ij <-> kl. Only the full matrix fflUj) is 
symmetric. Recalling the formula (no sum over indices) det(a l Mjjfe J ) = (JT,. a ibj) detM^, the 
determinant of (fTTj) is 

Here det C is a polynomial of order (n + l) 2 in the cosines, and it can be shown by direct 
computation with an algebraic manipulator that 

detC= ^Q G *^ detG ' ( 15 ) 

where Gij = cos 6^- is the Gram matrix with the convention Gu = —1, and G*y its cofactor 
matrix. Analogously one can show that all the minors are never singular and thus the rank is 
w( " 2 +1 ' > — 1 . This result shows that the non invertibility of Q has its root in the vanishing of the 
determinant of the Gram matrix on flat spacetime, namely in the closure of a flat n-simplex0 



2 On the other hand, the determinant of the Gram matrix for a spherical or a hyperbolic n-simplex does not 
vanish. This means that a perturbative expansion around these two types of background would directly yield 
an invertible quadratic order. We will comment on this point in the conclusions. 
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The elegant expression (TTTT) is fundamental to understand the connection between classical 
Regge calculus and the free propagator of the quantum theory: being the quadratic order of 
the expansion of the Regge action around a fixed background, the matrix Q is the key input 
to compute the free graviton propagator in Regge calculus, as we show below. In the special 
equilateral configuration, (TTTT) reproduces the results derived in [5] and it is used to study the 
free graviton propagator from spinfoams quantum gravity. 

The explicit expression flTTl) allows us to address two important questions on the hessian 
matrix: 

• What is its scaling? Recall that the free propagator scales as the inverse squared distance 
between the two points, evaluated with respect to the background geometry. In the 4d 
spinfoam calculations, these points are associated with centers of the triangles in the 
boundary of the 4-simplex. If the hessian Qijki = Qtt' is the crucial dynamical input as 
argued in [6], then it should encode the distance between the (centers of the) triangles 
t and t! . In particular, one might naively expect Q tt i to scale as some power of this 
distanced The general expression (TTTT) shows clearly that this naive expectation fails: 
even if Q has a simple geometry dependence, the various components of ( TTOT) can not be 
related to distances between points (such as centers of triangles) on the boundary of the 
4-simplex. 

Nevertheless we will see below that the matrix is indeed related to the free propagator: 
the right scaling only emerges when many 4-simplex contributions are glued together. In 
the setting of [6], this means that the contribution from the boundary state is crucial, as 
we discuss below. 

• What is the meaning of its non invertibility? We know that in the continuum the 
quadratic term of the Einstein-Hilbert action (to which the Regge action reduces in 
the continuum limit) is not invertible due to the invariance under diffeomorphisms. Can 
the non invertibility of Q have a similar origin? The null vector Vki does not appear to 
be related to a local symmetry such as the diffeomorphisms, but we will see below that 
it is indeed the remnant of the continuous diffeomorphisms on the single simplex. 

Notice that a natural way of making Qijki invertible is to add the rank 1 matrix built out 
of its null vector V^, 

~ 3(2-n) 

Qijki = Qijki + £V n VijVkl, (16) 

3(2-n) 

where the factor V « is chosen on dimensional grounds, and £ is a real number. 

In the next section we use a 3d model to show explicitly that Q is related to the free graviton 
propagator and that its non invertibility comes from the invariance under diffeomorphism of 
the original continuum action. In section [6] below we come back to the 4d case. 

3 Notice that this naive intuition could not be checked in the calculations of [6], because there only the simple 
equilateral configuration of the 4-simplex was considered, thus all the distances are the same and the scaling is 
trivial. 
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3 3d lattice model 



In this section we describe how the Regge action is related to the lattice graviton propagator 
developping the original analysis by Hamber and Williams [8]. We consider the 3d case where 
the lengths are the fundamental variables, thus we can address the questions at the end of 
the previous section without the extra complicacy of the constraints (jl]) needed in 4d. [We 
expect the results of this section to apply also to 4d length Regge calculus, and so to area 
Regge calculus if this is proved equivalent to length Regge calculus.] Using (TTTT) . we have for 
the tetrahedron 



1 i 



ijkl 



kl 



0V h 



ki 



9 sinfl,-,- V 2 



_ 



(17) 



with C the same matrix as before, and the following notation has been introduced: kl is the 
edge opposite to kl, and AjAj the (unique) pair of triangles different from A k A\. Notice that 
the need of the opposite edge in the columns of C comes from the fact that in (fTTj) we are 
deriving with respect to Vm = The double index notation is required only in this formula. 
For simplicity in the following calculations, let us introduce the following single index notation: 
we use l e = kj = Vjj for the edge e between % and j, and 6 e = 9jj for the dihedral angle associated 
with the edge e and thus hinging the pair of triangles Aj Aj. Finally, we define the tetrahedral 
matrix Q T ee ' = — fft = Q T TJki- This definition helps the comparison with the results of [H]. 

The 3d Schlafli identity reads ^2 e l e d9 e = 0, and the null vector is l e . As we discussed 
above for any dimension, the matrix (fTT|) does not have simple scaling properties with respect 
to the background geometry^ This matrix is the quadratic order of the Regge action on a single 

tetrahedron. To study how Q T is related to the free graviton 
propagator, we consider an infinite rectangular lattice, and 
divide each block into six tetrahedra as shown in the picture, 
namely drawing the diagonals on the six exterior faces plus 
the interior diagonal connecting the vertices and 7. All 
tetrahedra contain the vertices and 7, and can be obtained 
from the six possible ways of going from to 7 along three 
cartesian edges. Consequently each tetrahedron has the same 
edge structure: three cartesian edges, two face diagonals, and 

the body diagonal (07). All the edges are oriented (from smaller to bigger number), so that in 
the full lattice each edge can be identified giving its starting vertex and its direction, l e = l\, 
where e = (01) . . . (07) and v are the vertices in the lattice. Adapting cartesian coordinates 




1,0,0), (02) = r = (0,1,0), (03) 



(0,0,1), 



to this lattice, we have (01) = 
(04) = x^ + yv = (1, 1, 0), and so on. 

We now introduce the background around which we will compute the correlations. To 



4 The only elements for which we were able to find a simple relation to the distances are the ones for opposite 
edges: using ^ we have 



1 I 



9sin% V 2 



AiAj .2/1 

sin uji = IhItt 

3 6V 3 13 d 



1 



te),(y") sm< %Uy") 



where d 



(*j)>(y) 



is the distance between the midpoints of the two opposite edges Zy and hj, and 



angle between them. 



'mm 



(18) 
the 
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keep the lattice regular but also the tetrahedral matrix Q T as general as possible, we consider a 
rectangular background, with lattice spacings in the three cartesian directions given respectively 
by l"l = a, Z| = b and l\ = c, for any v. With this background, the six tetrahedra in each block 
are obtained by permutation of a, b, c, theyhave all the same volume abc/6, but the six edges in 
each tetrahedron are all different and thus the tetrahedral geometry rather general. Consider for 
instance the tetrahedron f = 0147: its edge lengths are (a, \J a 2 + b 2 , c, b, \/b 2 + c 2 , Va 2 + b 2 + c 2 ) 
and its Gram matrix is 



/ 



V 



-1 







-1 





ac 







Va 2 +b 2 Vb 2 +c 2 





-1 








Va 2 +b 2 



Va 2 +b 2 
-1 



where the numbering of row and column is 0,1,4,7 We use the following notation for the 
perturbative expansion: l e for the background values, £ e for the variables, and define the fluc- 



tuations <L 



l e , for which we also use the notation 81 where e denote an oriented edge. 



We are interested in computing the correlators between fluctuations around this background to 
quadratic order. To do so, we evaluate ( ITTj) on the background, for each of the six tetrahedra 
in the block. For instance, the tetrahedron f = 0147 gives 



1 

abc 



( 



\ 



o 

-b 2 (a 2 +b 2 ) 
b 2 (a 2 +b 2 +c 2 ) 

b 2 (a 2 +b 2 ) 
-{a 2 +b 2 )(b 2 +c 2 ) 



b 2 (a 2 +b 2 + c 2 ) 
-b 2 (a 2 +b 2 +c 2 ) 
-b 2 (a 2 +b 2 +c 2 ) 

b 2 {a 2 +b 2 + c 2 ) 






b 2 (a 2 +b 2 ) 
-b 2 (a 2 +b 2 +c 2 ) 
-fe 4 
b 2 {b 2 +c 2 ) 




-(a 2 + b 2 )(b 2 +c 2 ) 
b 2 (a 2 +b 2 +c 2 ) 
b 2 (b 2 +c 2 ) 
-b 2 (b 2 - 




-c 2 ) 



\ 



J 



where the numbering of rows and columns is (01) (04) (07) (14) (17) (47). Due to the regularity 
of the block, the matrices associated to the other tetrahedra can be obtained by appropriate 
permutations. Then adding up the matrices for the six tetrahedra one obtains the quadratic 
term of the Regge action for the block, 



- 8 e Q T ee >5 e i — — ^ 8 e Q h ee >8 e >. 



(19) 



r£b e,e'ST 



e.e'eb 



The equation above defines a matrix Q h associated with the block. This matrix is (19 x 19), 
corresponding to the 12 cartesian edges, the 6 face diagonals and the body diagonal. We 
report its explicit expression in the Appendix (see (1541) ). The total lattice action Sl is obtained 
summing over the blocks of the lattice, and it can conveniently written as a sum over vertices, 



v,v' e£n e'£v' 



(20) 



where e G v denotes one of the seven edges starting in v. For an infinite lattice, Q vv is an 
infinite dimensional matrix non vanishing only for first-neighbour vertices. Remarkably, we find 
for this matrix the following simple scaling, 



abc M * & ' 



(21) 
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where the entries of M are adimensional and independent on a, b, c; they are non zero only for 
first-neighbour vertices, and furthermore there are no correlations between edges starting at v 
and edges starting at v + x + y + z. The non zero entries of M vv are structured in 7x7 blocks 
reported in the Appendix. 

This simple scaling looks surprising (specially if compared with the non transparent scaling 
behaviour of the single tetrahedron matrix Q T ), and the reader might wonder if it is a pecu- 
liarity of the regular lattice chosen that would be lost on a non regular lattice. However, we 
show in Appendix [B] that the action (120]) can be obtained discretizing the quadratic term of 
the expansion of the Einstein-Hilbert action around the flat metric g^f = diag(a 2 , 5 2 , c 2 ) (see 
Appendix [B]). Remarkably, this discretization will produce exactly the scaling of ( l20i) . This is 
a peculiarity of 3d GR that we comment on in the Appendix. 

In spite of the appealing scaling, at this stage it is still unclear how this quadratic action 
will produce the correct scaling of the propagator. On the other hand, in the equilateral 
configuration a = b = c this block can be immediately recognized as the building block of 
the action term in [81, which is there shown to lead to the correct free propagator. In the 
following, we simplify^ and generalize the procedure of [5] to obtain the free propagator on this 
rectangular lattice. 

As in the continuum, the quadratic action (1201) has local symmetries that need to be gauge 
fixed in order to compute the free propagator. Calling Sg the gauge fixing term, we have 

W ee , = fUdSi 6 e 5 e , e-W-^lH (22) 

J i 

The standard procedure at this point would be to choose a particular gauge fixing, and going 
to Fourier variables to diagonalize the kinetic term. However, the use of a simplicial lattice 
leads to the presence of additional spurious variables in the action on top of the gauge degrees 
of freedom. This comes from the fact that the simplicial lattice gives seven directions for each 
vertex in the lattice, as opposed to the six variables per point of the metricj§ To isolate the 
spurious quantities, we make the following change of variables, 

^=E T ^ ; ' ( 23 ) 

0',e') 

where T vv is defined in the Appendix (see (1581) ) and has entries of length dimension one. The 
entries vanish for non first-neighbour vertices, thus this change of variables is local in the sense 
that it mixes only b\ which are first neighbours. The new variables h\ have length dimension 
2, and can be obtained discretizing the continuum metric fluctuations h^ u (x) with ihw and 
(fiu) i — > e according to the index correspondence (11) (22) (33) (12) (13) (23) i— > 123456. This 
index correspondence defines a matrix A^ u such that h\ = A^h^v), and selects h% as the 
spurious variable. In these new variables the action fl20l) takes the form 

sm = \ J2 ^^'^-^Ew) 2 . ( 24 ) 

v,v' V 
l<e,e'<6 

5 In particular, we do not use trace reversed variables as in [8]. 

6 We stress that this is an issue of the simplicial lattice alone and not of the theory being considered here, 
and indeed it would be present even if one was studying lattice gauge theory on a simplicial lattice. 



9 



where we see that hj completely decouples from the others. We can integrate out this spurious 
variable and ignore it in the following, that is from now on e = (1, ...,6)0 The explicit 
expression of L vv , which is rather lengthy, can be straightforwardly worked out from Q vv and 
T m given in the Appendix. Remarkably but not surprisingly, it can be also obtained by a 
discretization of the continuum Einstein-Hilbert action ^JgR around the background metric 
g hg = diag(a 2 , 5 2 , c 2 ), that is 




Ke,e'<6 



where hi = A^h^iv). The lattice derivatives are defined in the conventional way; denoting fi 
the unit vector in the direction /z and = (a, b, c) its length, we have A At 0(t>) = [0(t> + p) — 

L vv has six variables and three zero modes per vertex, that we also report in the Appendix 
(see fl62l) ). These three zero modes correspond to translations in the three cartesian directions, 
and are the remnants of the continuous diffeomorphism invariance on the lattice, as we discuss 
in the next section. This gauge symmetry can be fixed by adding a suitable term to f|2^|) . 
replacing the matrix L||,' by L~|/' = L||/' + L gi lf,. As L||, comes from the discretization of the 
continuum action, it is useful to construct also the extra term L gf ~|/ from the discretization of 
a continuum gauge fixing term. In particular it is convenient to choose the harmonic gauge as 
in [8], because it leads to a simple expression for the propagator @ The harmonic gauge fixing 
term is |\/ g hg C^C^ with = d p h fip — |(9 /J /z p p , where the indices are raised and lowered with 
the background metric. We report the discretization of C M in the Appendix (see (|63|) ). Using 
this gauge fixing, the matrix elements of L are given by 

Ltf,=abcY,D if Ay'D H , (25) 



where 
and 



D = diag^cT 2 , b~ 2 , c~ 2 , (a&)~\ (ac)^ 1 , (&c) _1 j , (26) 




-1 

-1 

10 

4 

4 

4 



A = ^ — + 5 v > iV -p, — 25 v i^ . (27) 



The quantity A can be immediately recognized as the laplace operator on the rectangular lattice 
here considered. Using the same [iv i— > e correspondence as above, it is easy to show that the 
kinetic term (1251) reduces in the continuum limit to the quadratic term in the harmonic gauge, 

StonM = l -jd i x^ V P " Upa A V> P^ P . = \ 9% 9 h ,l - \ 9ll 9% (28) 



7 This means that we do not look at correlations including the spurious variables. However, notice that such 
correlations vanish in the continuum limit (a, 6, c) h- > homogeneously. 

8 An interesting question, which we do not address here, is to consider other gauges such as the radial gauge 
suggested in [2T], and see how they affect the final expression of the lattice propagator. 
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(here A is the continuum laplace operator). 

At this point one evaluates the correlations in the new variables, 

/i bg „bg _|_ bg „bg _ o bg „bg 

obtaining the correct scaling with the distance d(v±, vn) between the two points and the correct 
tensorial structure of the 3d euclidean free propagator]^ We conclude that we recover the correct 
lattice propagator, generalizing the results of [8] to the rectangular lattice. The calculations of 
this section show how the single tetrahedron matrix (ITT)) adds up correctly to the lattice kinetic 
term of linearized Regge calculus. 

4 On the gauge fixing 

This connection with linearized Regge calculus allows us also to clarify the origin of the non 
invertibility of fflUj) . Let us go back to the three zero modes per vertex of the full lattice term 
L||, . These zero modes reported in the Appendix in ( l62|) are the three independent translations 
that can be performed at each vertex of the lattice, as shown also in [8] for the cubic case. This 
is the remnant of the continuous diffeomorphism invariance on a regular lattice, and consistently 
in the previous section we gauge fixed this symmetry simply discretizing the harmonic gauge 
for the continuum (linearized) diffeomorphism invariance. 

As an intermediate step towards the symmetries of the single tetrahedron matrix, we can 
look at a finite lattice built out of eight blocks. Such a lattice has an interior vertex. The 
quadratic action for this lattice (built as in (fl9l) adding up the tetrahedral matrices in the 
lattice) has the zero modes corresponding to the possible translations^ plus three additional 
zero modes corresponding to independent global rescaling of the lattice in the three cartesian 
directions. The latter are 

a 2 b 2 

4^4 + Nx a, 8y !->• 5y + Ny b, 5%+$ h-> 5 i+y + N & + N y (30) 

Va + b l Va + o 

and so on. Here N%, Ny- and N% are three free parameters. This symmetry is a consequence of 
the finiteness of the lattice. 

Next, we reduce the lattice to a single block, and look at the symmetries of the 19 x 19 
matrix (Tl9l . We find that all the zero modes related to translations are gone (the edge variables 
cannot be translated any longer), but the global rescalings (130]) are still present. These are 
the three zero modes of ffl9]) reported in the Appendix (see f]55]) ). Finally, the special case 
N$ = Ny = N% = iV gives the only zero mode that survives the initial tetrahedral matrix Q T , 
corresponding to the transformation 8 e i— > S e + Nl e . 

The analysis above allows us to make contact between the symmetry of Q T and diffeomor- 
phisms. It shows that the zero mode of Q T can be gauge fixed discretizing the continuum gauge 

9 We remind the reader the 3d GR without matter has no local degrees of freedom. Consequently this 
propagator is a pure gauge quantity, and does not represent a physical particle. 

10 Namely three zero modes corresponding to translations of the (edge variables starting at the) interior vertex, 
two zero modes for the 2d translations of the vertices on the exterior faces of the lattice, one zero mode for the 
Id translation of each vertex on the exterior edges of the lattice. 



11 



fixing term for the diffeomorphism invariance. For instance, one can use the harmonic gauge 
fixing for the full lattice described in the previous section to induce a gauge fixing term for 
Q T . Notice that this would produce a rather different term than the gauge fixing discussed in 
section [2] which, adapting ffTtj]) to 3d, reads 

ot = i^- (si) 

Indeed, the process of gauge fixing Q T alone (for instance using (I3"T|) ) and then adding up the 
gauge fixed matrices as in fl20|) to obtain a gauge fixed total action is ambiguous. This is 
because gluing tetrahedra together produces additional gauge degrees of freedom, of which a 
simple collection of copies of (|3T|) would not take care. This is evident from the presence in the 
gauge fixing term C (see the Appendix, (1631) ) of terms connecting edges belonging to different 
tetrahedra. From this perspective the use of gauges without derivatives, such as the radial 
gauge suggested in [21] . might be more promising when trying to relate the tetrahedral gauge 
fixing to the full lattice one] 11 ! 



5 General boundary formalism 

In the previous sections, we used a 3d model to show how the matrix Q T is related to the free 
graviton propagator, and how its non invertibility comes from diffeomorphism invariance. To 
make contact with the spinfoam calculations of the graviton propagator [6] , we use the general 
boundary formalism [221 ES]. As explained in [23] (and here adapted to 3d), we introduce a 2d 
closed surface E and split the integration of (T2"2"j) in the three regions of spacetime defined as 
the exterior of the surface, the surface itself, and the interior. In doing so we write the action 
as Sr = 5s + Scut where both pieces have a boundary term as described in (j2j. In particular 
we can choose E to be the surface of a single tetrahedron r, so that Sr, = S T + S out . We can 
then write (122]) as 

W ee , = fl[d8 t 5 e 5 e , e~ s ^ * T [5], (32) 

where i G r means the edges belonging to r, and the integration over the edges i ^ r outside 
r defines the boundary state 

* T [«J] = f Y[d6ie- S °^ 5] . (33) 

At this point we can proceed evaluating f]33l first, then plugging it into f]32l and evaluating 
the remaining integrals. As all we did is a splitting of the initial integration (]22l . we expect 
to obtain the same result of section [3] This formulation is nonetheless interesting because it 
seems more suitable for background independent approaches. In particular in [B] it is argued 

11 A way to characterize the lattice gauge fixing in terms of the null vectors, as it is done in ([31]) . is the 
following. If we denote by W the matrix constructed from the three null vectors (|62|) corresponding to the 
translations, one can show that \J g hg WC^C^W = A 1. This connection between the gauge fixing term and 
the null vectors could be used to identify the gauge fixing term in a sub block of the lattice. It suggests that the 
gauge fixing is dual to the null vectors, and so to replace (f3"Tj) by (,\/Vl~ l l~?~ . However the argument can not 
really be completed, as the discrete laplace operator appearing above is not well defined on a single tetrahedron. 
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that a quantity like (132]) - and its 4d equivalent for area variables, see below - emerges as the 
leading order of the spinfoam graviton propag ator@ As (|32|) gives the right propagator, this 
emergence would support the correctness of the semiclassical limit of spinfoams. However in 
[S] the boundary state \I/ T is not evaluated using the definition (1531) . but an ansatz for it is 
made. To shed light on \1/ T and improve the understanding of the ansatz and of the spinfoam 
calculations in general, we focus on (I3"3"]) in the rest of this section. 

The boundary state is a function of the fluctuations 5 e on the boundary of r. For the 
linearized theory, it is a gaussian in the fluctuations, and it is unique once the background and 
the gauge are fixed. It can be evaluated following the same procedure of the previous sections, 
with a few caveats that we stress below. The first difference concerns the expansion of S on t 
around the background. In expanding ( 122]) . we used the fact that the linear term = e e 
vanishes on the background solution, so that the leading order is quadratic. Yet this time, 
because both S T and S on t have a boundary contribution, the linear term does not vanish: it 
gives the extrinsic geometry of the boundary, according to For instance for S T we have 



where the dihedral angles 6 e of r give its extrinsic geometry, and the quadratic term is the 
matrix (flQl) . The expansion is analogous for S OVLt , in particular the linear term will be exactly 
the same (up to orientation), being r also the boundary of the outside region. At first sight 
the presence of linear terms might look like a big difference with the calculation of the previous 
section, but this difference is artificially created by the splitting of the original integration of 
(122]) in two parts as in (1321 . As (the background geometry of) r is adapted within a flat 
triangulation, the dihedral angles of these linear terms in fl33|) add up to give zero deficit angles 
after the integration over the lengths, and finally restore the structure of (122]) . This shows that 
it is natural to expect a linear term in the boundary state. 

A second caveat in the evaluation of the boundary state concerns how to deal with the 
spurious variables. In section 3 we isolated them performing the transformation (123}) . This 
transformation links only first-neighbours in the rectangular lattice, thus it is local on a box, 
but it is not well-defined on a single tetrahedron. This is not a problem as a tetrahedron has 
only six variables and so no spurious ones, but for a generic X one has to find a procedure 
to deal with this issue. In the special case when S is an infinite plane, the spurious variables 
can be isolated again by a change of variables such as (123]) . This case was considered in [24], 
where it was shown that the quadratic order of the boundary state for 4d length Regge calculus 
reduces in the continuum limit to the correct expression (given for instance in [251 26J ) . 

The last caveat concerns the gauge fixing. Evaluating the integrals defining (13"3"|) requires a 
gauge fixing, thus the boundary state is gauge dependent. This point is particularly relevant 
if one is given the explicit boundary state (for instance through an ansatz) and plugs it into 
( 132]) to compute the propagator: to get the correct answer, the remaining integrations in ( |32|) 
have to be gauge fixed consistently with the choice of gauge for the boundary state. Notice 
that contrarily to the continuum case, the gauge chosen might not be transparent in the lattice 
explicit boundary state. The most natural procedure to avoid this problem is to split the full 
lattice gauge fixing term into boundary and bulk contributions. 

12 In particular, it emerges from the simplest spinfoam, the one characterized by a single vertex. More general 
spinfoams can be put in correspondence with different choices of E. 




(34) 
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6 Conclusions and perspectives on the spinfoam 
calculations of the graviton 

In this paper, we studied the connection between the hessian matrix of the Regge action for 
an n-simplex and the free propagator of the linearized quantum theory. Our first result is the 
general formula for the hessian matrix, given in (ITTj) . This is needed to compute the perturbative 
expansion of the Regge action around an arbitrary background, and enters the evaluation of 
the free propagator, or 2-point function. The latter scales as the background distance between 
the two points to the power 2 — n, so one might naively expect a similar scaling for the hessian 
matrix. Contrarily to these expectations, the hessian shows no transparent scaling (with the 
only exception described in footnote H]). Furthermore, we showed that the hessian has a single 
zero mode, and thus is not invertible. We computed explicitly the determinant of the hessian, 
equation (fTST) . and showed that the zero mode has its root in the flatness of the n-simplex. 

To understand both the scaling and the non invertibility of the hessian, we considered a 
3d lattice model. We constructed the full lattice action adapting (ITTj) to a rectangular lattice. 
Then we generalized (and simplified) the techniques introduced in [7J and [8] and show that 
upon properly gauge fixing the action one obtains the correct free lattice graviton propagator, 
given in ( |29l) . The procedure answers the issues raised above in the following way. First, the 
correct distance scaling, hidden in the single tetrahedron term Q T , emerges when the full lattice 
is considered (from the point of view of the general boundary proposal ( 1321) . this means only 
when the correct boundary state is added). Second, the zero mode of the hessian is explicitly 
related to the diffeomorphism invariance of the continuum action. In particular, we showed 
that the hessian matrix can be gauge fixed discretizing the continuum gauge fixing term for 
the diffeomorphisms. A consistent gauge fixing is crucial to obtain the correct free graviton 
propagator. 

These results were obtained explicitly in 3d, but can be extended to 4d length Regge calculus. 
On a hypercubical lattice, it was shown in [7J that linearized quantum Regge calculus gives the 
correct free spin-2 propagator in the continuum limit. To relate those results to the hessian, 
one can use the 4d Q a to construct the full quadratic action for a 4d rectangular lattice, where 
the 4-simplices will have all different edge lengths. Then one can proceed as in section [3] to 
generalize the results of [7j to this rectangular lattice. The calculations become lengthier, but 
there are no technical obstructions. 

The situation is different if one uses the areas as fundamental variables, as suggested by the 
spinfoam calculations. Studying linearized area Regge calculus (and its propagator) is a more 
complicate problem, due to the lack of knowledge about the constraints needed to treat the 
areas as independent variables, as discussed in section [2J An intermediate step considered in 
[HI ED] is to take the functions of the edges, A t (l e ), namely one assumes to have solved 

the constraints. Accordingly one looks at dependent fluctuations 8t(l e )- The resulting theory 
can be shown to be equivalent to conventional length Regge calculus. In order to extend our 
conclusions to area Regge calculus, it is crucial to study the constraints in Due to the lack 
of knowledge about their explicit form, we are not in a position to draw for area Regge calculus 
the same conclusions we have for the conventional theory. 

Nevertheless the results here obtained help clarify important aspects of the spinfoam gravi- 
ton calculations of [6], in particular the issue of gauge fixing. Let us summarize the steps from 
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( |32|) to the final propagator expression (|29|) : (i) we expand S T to second order around the 
background; (ii) we gauge fix its symmetry; (iii) we evaluate the boundary state (I33I) . again at 
quadratic order, in the same gauge; (iv) the remaining integrations in (|32|) give the correct free 
propagator. 

If we compare these steps with the general setting of the calculations in [5], the step (i) 
is the same, but the steps (ii) and (iii) are not performed. Concerning (ii), there is no gauge 
fixing of the symmetry of S T . Concerning (iii), the boundary state is not evaluated, but rather 
a gaussian ansatz is made. This gaussian ansatz has a linear and a quadratic term. The 
linear term is perfectly consistent with the one found in the previous section, although in [6] 
its presence is justified on different groundsill The quadratic term, on the other hand, is 
left generic, with free parameters giving the spread of the gaussian. As a consequence, those 
calculations can not be used to check the correct distance dependence of the propagator, nor its 
tensorial structure, because both aspects rely heavily on the precise structure of the boundary 
state. Based upon a generic boundary state, the calculations can only test the correct overall 
scaling of the propagator, and this is the result indeed achieved in In principle, one can 
impose the correct result to constrain the free parameters entering the gaussian. However 
this is rather subtle for the presence of the gauge fixing term in the boundary state. In fact, 
concerning the absence of gauge fixing in [6], point (ii), the reader might wonder why those 
calculations achieve a finite answer. The reason is simply that the genericity of the boundary 
state acts as a gauge fixing term for Q T , making the overall (boundary plus bulk) quadratic 
term of the action invertible. Therefore the free parameters appearing in [6] should be fixed 
not only by a dynamical derivation of the ansatz, but also by making a precise choice of gauge 
for the full lattice. This is an unusual gauge fixing procedure: in order to get the right answer, 
the boundary state must carry also the correct additional gauge fixing information for the bulk 
term. A more natural procedure that emerges from the results presented here is to first choose 
a gauge for the boundary state, then evaluate it through its definition (1551) (or alternative 
ways 14 !) . and finally fix the symmetry of the bulk action in the same gauge. 

Since the gauge fixing is hidden in the boundary state it seems important to us to explore 
the dependence of the boundary state on the gauge fixing in a more systematic manner at least 
in Regge calculus. The setup which we presented in this work appears to be a good starting 
point for this task. Possible extensions which are relevant for the spinfoam calculations are the 
use of different gauge fixings and even less regular lattices. 

Understanding these issues for the graviton is particularly relevant since we are computing 
a gauge dependent observable. Another strategy to study the semiclassical limit of spinfoams 
could be to look at gauge invariant observables and compare them with their counterpart in 
Regge calculus. For instance it would be interesting to adapt to this discrete Regge setting the 
framework to study GR perturbation theory in a gauge invariant way developed in [27]. 

A delicate issue to be understood in the general context concerns the spurious variables. As 
we have seen identifying them is crucial to reproduce correctly the propagator. The spurious 

13 Notice that in [6] the linear term is imaginary, thus giving a phase to the boundary state. The factor i that 
makes the linear term into a phase can be understood recalling that in the spinfoam calculations one is using a 
kernel exp{— iS} as opposed to exp{— S} used here. 

14 It is useful to remind the reader that, at least in the conventional context of background dependent quantum 
field theory, the boundary state can be directly extracted from the kernel. This procedure for linearized gravity 
is described in [55] , 
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variables come from the simplicial nature of the lattice, and have nothing to do with gauge 
variables, thus fixing the gauge does not take care of them and they have to be identified in 
an independent way. In the regular lattice considered here we were able to do so studying the 
special way 8% enters the action, but on a general irregular lattice such as the one arising in a 
general spinfoam model this identification could be a very challenging issue. 

Another interesting point that could be further developed is the perturbative expansion 
around a constant-curvature background. As it is clear from ( [15]) . the hessian becomes invertible 
if the background is a spherical or hyperbolical n-simplex (respectively positive or negative 
Gram determinant). We leave this point open, and make the simple following remark. As 
we are dealing only with riemannian signatures, these backgrounds can be associated with the 
perturbative expansion in the presence of a non-zero cosmological constant A. In fact, A > 
(respectively A < 0) can be associated with a global spherical (hyperbolical) topology. Then 
this spacetime could be realized on a lattice made out of spherical (hyperbolical) n-simplices. 
In this case, the n-simplex hessian is invertible, and gauge degrees of freedom in the full lattice 
action only arise from the gluing of the n-simplices. 
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A Explicit evaluation of Q 

In this Appendix we report the calculations leading to (fTTT) . It is convenient to use the affine 
(or barycentric) coordinates introduced in [28] for Regge calculus (see also [29]). These are 
defined as follows. Let Vj, j — 1, . . . , n + 1 be vectors from some arbitrarily chosen point in M. n 
to the n + 1 vertices of an n-simplex a. Then 



1 d 

^^^-——^tk (35) 
it=i 

define the affine basis. Notice that Ylj e j = 0j so this basis is overcomplete. As a consequence, 
the affine coordinates x 3 of a vector x = 5^ ■ x 3 ej are not unique. Uniqueness can be obtained 
requiring J2j & = 0- 

We can also introduce a dual affine basis e 3 defined by 

ej ■ e k = S k = 5 k - -J— (36) 
3 3 3 n+ 1 

where 5 k is the Kronecker delta. The dual basis satisfies e J = 0. This enables us to define 
affine components of an arbitrary tensor in IR n by contracting it with the appropriate product 
of affine basis vectors ej, e k . 
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The key quantity for our calculations is the metric tensor in these coordinates. It can be 
shown that the affine components of the metric tensor are [28] 



9ij .) 



i£W- (37) 



In particular one can easily check that contracting (157|) with the edge vectors (e i:) ) fc = Sf — 6* 
gives indeed the length I 2 - of these vectors. Laplace's formula for the determinant of gives 
the squared of the n-volume of the simplex, 



1 



v -"T~m e e 9kih ' ■ ■ 9k n i n , [oo) 



n! 



where the affine epsilon tensor is defined as 

~j v .. jn _ | + 1 if the permutation j x ■ ■ -j n j n+ i is even 
1—1 if the permutation j\ • • • jnjn+i is odd 

and is vanishing if {ji • • • j n } can not be completed to a permutation of {1, . . . , n + 1}. 

The inverse affine metric components are defined by gug l% = g %l gm = 8 % k , and one can 
use ( 1581 to write 

= (^^^^^•••ffW. = ( 4 °) 

The diagonal components of the inverse metric g w are proportional to the (n — 1)- volumes of 
the (i)-subsimplices, 

1 V 2 

r = -ir 2 - (41) 

Analogously, the (n — 2)-volumes Vij of the (2j')-subsimplices (i.e. in 4 dimensions the areas of 
the triangles not containing the vertices i and j) can be obtained by fixing two indices in both 
epsilon tensors of the RHS of (jlD]L 

V 2 = - 7ijk3-k n ~ijh-l n ~ . . .n, , (AO) 

ij ((n-2)!) 3 33 
From this formula we can immediately write the derivatives with respect to the edge lengths, 



dV 2 (n-2 ( 



l± ~ ^ )_^3hk A ...k n ^mU-l n . . . ~ knin _ (43) 



«L ((^-2)!) 3 

The inverse affine metric can be used to obtain the matrix (fill of derivatives of the dihedral 
angles with respect to the edge lengths. To do so, we need to relate the dihedral angles to the 
inverse affine metric. The dihedral angle 8^ is defined to be the angle between the (inner) 
normals n(i), n(j) to the faces i and j. The affine components of these are given simply by 
n{i)k = $1, and their norm is 



2 = J2n(i) k n(i)rg kl = = t = (44) 



l kl 
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where in the last step we used (l4Tj) . This gives for the dihedral angles 

-s% = -^%#^ = -^£=- (45) 

Notice the minus sign, due to the fact that the 9 are defined as the inner dihedral angles. This 
is consistent with the convention Ga = — 1 for the Gram matrix used in section [3 Using (T4"TT) 
and (PIS]) , we can write the affine inverse metric as 

p = _±vmfi m9v _ (46) 



It is now straightforward to obtain the derivative of the dihedral angles with respect to the 
length variables. First, we take the derivative of equation (|45p . getting 



d9kl_ = Jkrn 1 = / khml kmhl _ El-khgkrn _ 9^~lh~lu\ (4?) 

dl hm sm9 M ^kT~u\ y y y g kk " y 9 U J 

Next, we use ( 1461) to replace the inverse affine metric components by functions of the dihedral 
angles, obtaining the final expression in ( fill) : 

ddki 1 Ihm VhV m , 

COS 9 kh COS m i + COS flfc™ COS V M + 



9/ ftm n 2 sin0 w V 2 

cos9 k i(cos9 kh cos9 km + cos 9i h cos 9i m )^. (48) 

Note that this formula holds for general dimensions n. 

We also describe an alternative derivation of the formula, which gives the intermediate 
step in ( ITT]) , and allows a more geometrical understanding of the different contributions to the 
differential of the dihedral angles. We start from the general law of sines in any dimension (a 
derivation can be found at the end of this section), 

n V-V 

sin 9a = / 1J . (49) 



Taking the total derivative of equation (1491) we obtain 

n 1 VVij fdV dVij dV t dVj 
~ 7Z. TT ^ZTE~ 7777" + 



(n-l) cos% ViVj V. K Vij Vi Vj 



= \ tan% (l m + P{ijf m - P{if m - P{ 3 ) k m ) ~g ml dg H - (50) 

In the last line we used that the volume differentials of the simplex and the subsimplices can 
be written as 

^ = \H^ k m 9 ml dg kl (51) 
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where x can stand for a (then V a = V) or i or ij. P(x)^ n for x = i and ij is the projector in affine 
coordinates onto the subsimplices (i) and (ij), respectively. For x = a we have P(a)^ n = 8^. 
These projectors can be expressed in terms of the inverse affine metric components, which leads 
again to equation (jUJ). From the last line of fl50|) one sees easily that the total derivative of the 
dihedral angle vanishes if contracted with a variation vector 8 corresponding to a global scaling 
of the simplex, i.e. such that dcjki(8) oc g ki . 



The law of sine (l4"9l can be derived by the following consideration [3U]. Let d(i) be the 
distance from the vertex i to the subsimplex (i) and d(i\j) be the distance in the subsimplex 
(j) of the vertex i to the subsimplex (ij). Then we have a right triangle where two of the sides 
are d(i) and d(i\j) and the angle opposite to the side of length d(i) is the dihedral angle %. 
Hence we have 



sin 8ij 



d(i) 



n 



VVi 



ij 



ViVj 



d(i\j) (n-l) 

where in the second step we used that for any m-simplex a m we have 



(52) 



d(i\a n 



m\ V(a r 



(m - 1)! V(i\a m ) 



(53) 



Here V(a m ) is the volume of the simplex a m and V(i\a m ) is the volume of the subsimplex of 
a m , that one obtains by removing the vertex i from a m . 



B 3d lattice model 



In this Appendix we report explicitly the matrices used in section [3] and give more details on 
the calculations. We begin reporting the 19 x 19 block matrix defined in (f!9l) . Adding up all 
the tetrahedral contributions, we find 



Q 



abc 



(54) 



where 
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/ 



The ordering of rows and columns is (01, 02, ... , 07, 14, 15, 17, 24, 26, 27, 35, 36, 37, 47, 57, 67), 
where ij denotes the edge between the vertices i,j £ (0, 1, . . . , 7) of the block. We find that 
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Q h ee' has three null vectors v\, i = x,y, z, given by 

hv% = (1,0,0,1,1,0,1,0,0,0,1,0,1,1,0,1,0,0,1), 
hv\ = (0,1,0,1,0,1,1,1,0,1,0,0,0,0,1,1,0,1,0), 
hvl = (0,0,1,0,1,1,1,0,1,1,0,1,1,0,0,0,1,0,0). (55) 

These null vectors where discussed in section HI and correspond to the scaling symmetry fl30|) . 

The next step in ( |20l) is to rewrite the total action as a sum over vertices. For each pair of 
vertices {v,v') we get a 7 x 7 matrix, with the same scaling as ( 15^1) . encoding the correlations 
between the seven edges starting from the vertex v and the seven edges starting from the vertex 
v' . We write this as in (|2~T|) . 

and report here the relevant elements of M~~, . Denoting 5 v i >v the Kronecker delta for the lattice 
vertices, we have 

Mvv' q r A/f vv ' A j_ A ft/T vv ' A _|_ A l\/r vv ' A — A 

^ — £U v > jV7 lvl l2 — u v' ,v+x+z T Vv',v—y—zi lvl 13 — u v' ,v+x+y t Uv\v— y — z , lvl \4 — u v' ,v u v / 1 v—y^ 

A^15 — &v',v 3v',v — z, -^16 = &v',v+x Sv' ,v—y — z> ^17 = ^v',v — y ~\~ $v' ,v — zi 

Mvv' r ■ r ]\/T vv ' 9 A 1\/T VV ' A _i_ A l\ f vv> A A 

21 — Vv' ,v— x — z T" ^tj'jV+y+i, Jw 22 — zc, u',t?i Jvi 23 — °V ,v+x+y i Uv',v — x — £ , Jw 24 — °v' ,v 

-^25 = fiv\v— x— z: -^26 = ^v',v ^v',v— -^27 ~ 3v',v — x ^v',v — &•> 

Mvv' _ _ r r ]\/f vv ' — _A _ A , - A/7"" 1 '' — A - _ A . . H/T vv ' — A 

41 — "ti'.u Vv',v+yi lvl 42 — °«',t) "u.'u'+cc ) iu 43 — "t>',ii — z "u'^+x+i/! iw 44 — ^^v',vi 

— ^u',u— i 4" <^ti',ti+i/, -^46 = 3v',v — i H~ ^v'.v+Sj -^47 = ii 

-^71 = Sv' ,v+y 5-v' ,v+z, — <5'u',u+£ + fiv',v+zi ^73 = ^v' ,v+x ^>v' ,v+y i ^74 = _i ^d',i; ~ ^»',u+Z] 

-M75 = Sv',v 3v',v+p> ^76 = <^u',t> ^v',v+xi -^77 = 2(5^'^. 

The remaining non zero entries can be found following the pattern given above. 

To compare these results to the continuum Einstein-Hilbert action it is convenient to change 
variables from the 5 e labeling the fluctuations of the edge lengths to (suitable discretizations 
of) the metric perturbations = — g^f. The latter can be discretized using the lattice 
cartesian coordinates x M = (1,0,0), etc., introduced in section [3) The change of variables can 
be found expanding to first order in the 5 variables the following relations, 

g^x^ x v = (h + 5i) 2 

% v {x + yf (x + vY = {h + h) 2 

9^(x + y + zY{x + y + zY = (l 7 + 5 7 f . 

These relations parametrize the six metric components per point in terms of seven edge lengths 
per lattice vertex. This brings us immediately to the problem of the spurious variables: one of 
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the seven 5g is redundant. In principle we can define a degenerate transformation 5\ i— > h^(x), 
but it is more convenient to introduce an auxiliary seventh metric variable per vertex, h%. We 
can then introduce discrete metric variables /i| i— > (h^ix), h%), where for the first six variables 
we have the mapping jh« and (11) (22) (33) (12) (13) (23) i— > 123456. Because of the absence 
in (155]) of correlations between edges starting at v and edges starting at v + x + y + z (this can 
be seen from the lack of terms like 5 v r tV+ x+y+z in it is convenient to relate this spurious 

variable to the fluctuation Sy, requiring that 8\ . . . 8q are functions of the six metric variables 
only. This can be achieved with the following transformation, 



v',e' 



(57) 



with 



/ 



2aS v r !V 



-aS v ',v 
d5 v > ^ v 






2bS v / 


-bS v , 


-b<5„ 



2cS 







-cS v 
-cS„ 



(58) 



\ 





Va 2 + b 2 S v i >v 

Va 2 + c 2 5 v j iV 

Vb 2 + c 2 S v , tV 



-^Va 2 + b 2 (8 v ,, v + 8 V , <V+B ) -^Va 2 + c 2 (8 v , :V + S v i >v+f) ) -|Vb 2 + c 2 (S v \ v + V,u+a) Va 2 + b 2 + c 2 5„/,u J 

Notice that T vv> has length dimension one, and thus h\ length dimension two. This change of 
variables leads to a decoupling of the hy variables: indeed, in these new variables the action 
reads 

Sm = \ £ KL%h%-—jy^. (59) 

v,v f V 
l<e,e'<6 

Remarkably, it can be proved that the first summand gives the discretization of the second vari- 
ation of the Einstein-Hilbert action y/gR around the background metric g hg = diag(a 2 , b 2 , c 2 ), 

[V; g hg ) = v 7 ^ {-{h^drhr + Ih^d^br - \h^d v K + \h»d v d v K) . (eo) 

Furthermore, an explicit calculation shows that 

5 ( 2 )[V;5 bg ] = ^^ 2 )[V;^]- (6i) 

where rj = diag(l, 1, 1) is the euclidean metric. This somewhat surprising result shows that 
the simple scaling of (I2TI) (or (!56j) ) is also a property of the continuum action (the additional 
lengths scalings present in (121!) and (1561) are absorbed into the transformation (|57|) from the 
S- variables to the metric variables). We stress that this is peculiar to three dimensions: the 
same action ( IBTil) in 4d does not have the property (T6TT) . As we comment below, this 3d scaling 
( I6"TT) is related to the fact that three dimensional gravity has no local degrees of freedom. 
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The remaining step is the analysis of the gauge symmetries, and their fixing. In terms of 
the 5-variables, the total lattice kinetic term (|56|) (and so ([2Tp in the main text) displays three 
null vectors per vertex, given by 

\$v' ,v+x $v' ,vi 0; 0) &v' ,v+x+y $v' ,Vi &v' ,v+x+z $v' ,Vi 0; &v' ,v+x+y+z $v' ,v)i 
le (0, $v' ,v+y ^v',v^ 0? &v f ,v+x-\-y ^v f ,vj O5 ^v' ,v-\-y-\-z &v',vi &v' ^v+x+y+z 

lg Wg (0, 0, 8 v '^ v -\-z &v' ,Vi 0> &v' ,v+x+z $v' ,Vi &v' ,v+y+z $v' ,Vi &v' ,v+x+y+z ^v',v}- (62) 

These null vectors correspond to the induced change of edge length if a vertex is translated in 
either the x, y or z direction. Therefore (12T|) has a gauge symmetry given by translations in 
3d space, which is what remains of the continuum diffeomorphism symmetry on the regular 
lattice. 

Understanding these symmetries in terms of their continuum counterpart helps choosing 
a gauge in which the continuum limit of the graviton propagator looks more transparent. 
We know from the continuum calculations that the propagator is simplest in the harmonic 
gauge Cfj,(x) = d p h fip — \d^h p p = 0. This quantity can be immediately discretized using the 
correspondence with the h\ variables introduced above, and one obtains C M (x) t— > Yle^e^l 
where 

(2a^( &v' ,v &v',v-\-x) ^z(.$v' ,v $v' ,v+x) 2<j? i^v' ,v ^i>'.r— x) 
2^?(.^v',v ^v' ,v+z) 2fiz{,$v' ,v $v' ,v+z) 2c? ( ^v' ,v &v',v+z) 

~gl(3v',v &v',v—y) ~£l(3v',v $v',v— z) \ 

j?(fiv',v - $v',v-x) ^(S v ',v - $v',v-z) ■ (63) 

^($v',v ~ 3v',v-x) ]p(dv',v ~ $v',v-y) J 

Analogously, one can discretize the gauge fixing term | a/ g hs C^C^ leading to ( |25i) in the main 
text. This completes the derivation of the main results presented in section [3j 

To conclude, notice the non trivial scaling of (1631) . and in particular that it does not agree 
with the simple scaling of fl61~]) . This non trivial scaling is crucial in order to obtain the correct 
scaling of the discretized laplacian in (J2j 



A = (S v ',v+ii + 8v',v-fi - 2Sv',v)/ll- (64) 

The fact that the non trivial scaling of f)64p comes from the gauge fixing terrrf^l and not from 
the kinetic term of the action might look confusing at first. However, this is a well known feature 
of 3d general relativity: it is a theory with no local degrees of freedom, and the propagator we 



15 Indeed if one added a term C iL C v rf' v /abc which has the same scaling of (|61j) one would end up with the 
laplacian for the background metric rj = diag(l, 1,1), namely 

A = {pv,v+p, +8v',v-p, - 26„',„). (65) 
A 

Notice that of course C u r]^ v / abc is not covariant (it has the wrong background metric) thus it is not an 
admissable gauge fixing term. 
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are computing is a pure gauge quantity] 16 ! Consistently, its dynamical properties are completely 
determined by the gauge. In 4 dimensions on the other hand flSTl) is no longer valid, hence the 
above mechanism does not apply. 
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